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Abstract 



We consider the single server queue with service in random order. For a large class 
of heavy-tailed service time distributions, we determine the asymptotic behavior of the 
waiting time distribution. For the special case of Poisson arrivals and regularly varying 
service time distribution with index —v, it is shown that the waiting time distribution is 
also regularly varying, with index 1 — v, and the pre-factor is determined explicitly. 

Another contribution of the paper is the heavy-traffic analysis of the waiting time 
distribution in the M/G/l case. We consider not only the case of finite service time 
variance, but also the case of regularly varying service time distribution with infinite 
variance. 
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1 Introduction 

We consider a single server queue that operates under the Random Order of Service discipline 
(ROS; also SIRO = Service In Random Order): At the completion of a service, the server 
randomly takes one of the waiting customers into service. Research on the ROS discipline has 
a rich history, inspired by its natural occurrence in several problems in telecommunications. 
The M/M/l queue with ROS was studied by Palm [25J, Vaulot [28] and Pollaczek [261127]. 
Burke jTT] derived the waiting time distribution in the M/D/l case. An expression for 
the (Laplace-Stieltjes transform of the) waiting time distribution for the M/G/l case was 
obtained by Kingman [20] and Le Gall [22J; the former author also studied the heavy-traffic 
behavior of the waiting time distribution, when the service times have a finite variance. Quite 
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recently, Flatto |15] derived detailed tail asymptotics of the waiting time in the M/M/l case. 
As pointed out in Borst et al. [Bj, this immediately yields detailed tail asymptotics of the 
sojourn time in the M/M/l queue with Processor Sharing, because these two quantities are 
closely related in a single server queue with exponential service times. 

In the present study, we are also interested in waiting time tail asymptotics of single 
server queues with ROS. However, here we concentrate on the case of heavy-tailed service time 
distributions. The motivation for this study is twofold. Firstly, an abundance of measurement 
studies regarding traffic in communication networks like local area networks and the Internet 
has made it clear that such traffic often has heavy-tailed characteristics. It is therefore 
important to investigate the impact of such traffic on network performance and to determine 
whether possibly adverse effects can be overcome by employing particular traffic management 
schemes. One possibility is to modify the 'service discipline' (i.e., scheduling mechanism); this 
may lead to a significant change in performance |7J. 

Secondly, in real life there are many situations in which service is effectively given in 
random order. Our own interest in ROS was recently revived in a joint project with Philips 
Research concerning the performance analysis of cable access networks. Collision resolu- 
tion of user requests for access to the common transmission channel is being handled by a 
Capetanakis-Tsybakov-Mikhailov type tree protocol [3J. That collision resolution protocol 
handles the requests in an order that is quite close to ROS [9] . 

We now present an outline of the organization and main results of the paper. Section [2] 
contains preliminary results on the busy period and waiting time tail behavior in the GI/G/1 
queue with a non-preemptive and non-idling service discipline. They are used in Section [3] 
to study the waiting time tail for the GI/G/1 queue with service in random order. The tail 
of the service time distribution is assumed to be in the class Cf]T>. This class contains the 
class of regularly varying distributions; these two classes, and others, are briefly discussed in 
Appendix A. We sketch a probabilistic derivation of the asymptotic behavior of the waiting 
time distribution, deferring a detailed derivation to Appendix [Cj For large x, P(M /7 ros > x) 
is written as a sum of four terms, each of which has a probabilistic interpretation. These 
interpretations are based on the knowledge that, for sums of independent random variables 
with a subexponential distribution, the most likely way for the sum to be very large is that 
one of the summands is very large (similar ideas were developed in [2] for a class of stochastic 
networks - see the so-called 'Typical Event Theorem' there). For example, the first of the four 
terms equals p times the probability that a residual service time is larger than x, p denoting 
the traffic load. The probabilistic interpretation is that one possibility for the waiting time of 
a tagged customer to be larger than some large value x is, that the residual service time of the 
customer in service upon his arrival exceeds x. The other three terms are more complicated, 
taking into account possibilities like: A customer with a very large service time has already 
left when the tagged customer arrived, but it has left a very large number of customers behind 
- and the tagged customer has to wait for many of those (and newly arriving) customers. 

In the subsequent sections we restrict ourselves to the case of Poisson arrivals. In the 
case of an M/G/l queue with regularly varying service time distribution, we are able to 
obtain detailed tail asymptotics for the waiting time distribution, in two different ways: (i) 
in Section [4] we apply a powerful lemma of Bingham and Doney |4] for Laplace-Stieltjes 
transforms (LST) of regularly varying distributions to an expression of Le Gall [22] for the 
waiting time LST in the M/G/l queue with ROS; (ii) in Section [5] we work out the general 
tail asymptotics of Section [3] for this case. Either way, the waiting time tail is proven to 
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exhibit the following behavior in the regularly varying case: 



P(W BOS >x)~-B-h{p)P(Bf w >x) 1 x^oo. (1.1) 
1-p 



Here, and throughout the pa 



in Formulas (4.11) and (4.12 



iper, f(x) ~ g(x) denotes lim^oo f{x)/g(x) = 1; h(p) is specified 
Bf w denotes the forward recurrence time of the service times, 
i.e., the residual service time. It is well-known that, with B denoting an arbitrary service 
time, 

P(B^ > x) = r P{B > U) &U, x > 0. (1.2) 
Jx EB 

Note that, except for Poisson arrivals, B^ w has a different distribution than the residual 
service requirement of the customer in service at arrival epochs. 



Formula (1.1) should be compared with the waiting time tail asymptotics in the M/G/l 
FCFS case [21]: 

P(W / FC FS > x) ~ -^P{B fw > x), x^oo. (1.3) 
1-p 

We shall show that h(p) < 1, which implies that ROS yields a (slightly) lighter tail than 
FCFS. 

In Section [6] we allow the service time distribution to be completely general. We study the 
waiting time distribution in the case of heavy traffic (traffic load p "f 1). When the service time 
variance is finite, we retrieve a result of Kingman [20J. When the service time distribution is 
regularly varying with infinite variance, we exploit a result of [S] to derive a new heavy-traffic 
limit theorem. 

The paper ends with four appendices. Appendix |A| discusses several classes of heavy-tailed 
distributions. Appendices [B] and [C] provide the proofs of two theorems. In Appendix [D] we 
state and prove a lemma that is not explicitly used in the paper. However, it has been very 
useful in guiding us to the proofs of our main results. Essentially, the lemma states that when 
interested in events involving a large service time, we may in fact ignore the randomness in 
the arrival process. 

Remark 1.1. A different way of randomly choosing a customer for service is the following. Put 
an arriving customer, who finds n waiting customers, with probability ^j-j- in one of the posi- 
tions 1, 2, . . . , n + 1, and serve customers according to their order in the queue. Fuhrmann and 
Iliadis [T7] prove that this discipline gives rise to exactly the same waiting time distribution 
as ROS. 



2 Preliminaries: Busy period and waiting time 

We first focus on the busy period of the GI/G/1 queue. For the time being we may take the 
service discipline to be the familiar FCFS, since the busy period is the same for any non-idling 



discipline. At the end of this section - in Corollary 2.5 - we use the results on the busy period 
to state a useful relation for the waiting time in any non-idling service discipline. 

Let us introduce some notation. The mean inter-arrival time is denoted with a and the 
random variable B stands for a generic service time, with mean EB = 0. A generic busy 
period is denoted with the random variable Z and r is the number of customers served in a 
busy period. The residual busy period as seen by an arriving customer (i.e., the Palm version 
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associated with arrivals) is denoted with Z rp . As before, we use B^ w to denote a random 
variable with the forward recurrence time distribution of the service times. 

For the GI/G/1 queue the proof of the following proposition is given in [16J. The defini- 
tions of the widely used classes S* , ITZV, C and T> can be found in Appendix |A} The first 
proposition can be specialized to the M/G/1 queue by substituting Er = 

Proposition 2.1. If B £ S* , then, for any < c\ < 1 < 02, 

P(r > n) 

lim SU P C . D/D ^ I „Jn IvT < 1, (2-1) 

n- 

and 



P(t > n) 


EtP(B > ana(l 


-P)) 


P(t > n) 




EtP(B > c 2 na{\ 


-P)) 


P(Z > x) 




EtP{B > d ia ?(l - 


~P)) 


P(Z > x) 




EtP(B > d 2 x(l - 


~P)) 



lim inf - — y ' rr > 1. (2.2) 

Similarly, for any < d\ < 1 < d 2 , 

lim SU P t, D / d Z\ ^ x i ( 2 - 3 ) 

X 

and 

lim inf — ' v " " 7 ^ > 1. (2.4) 

In particular, if B G T7^.V i/ien 

P(r > n) ~ ErP(S > na(l - p)), as n -)• oo, (2.5) 

and 

P(Z > x) ~ ErP(B > x(l -p)), as x -»■ oo. (2.6) 

The next proposition gives the asymptotics of the distribution of the residual busy period. 
Heuristically speaking, it indicates that a large residual busy period requires exactly one 
large service requirement (in the past). When analyzing waiting times (and residual service 
requirements) this result proves to be very useful as we shall see later. In fact, we shall 



sharpen the statement of the proposition (in line with the heuristics) in Corollary 2.3 
Proposition 2.2. If B £ £f)T>, then 

oo 

P(Z r P>x) ~ Yl P i B ~m > {x + ma){l - p)) (2.7) 



m=l 
P 

1-p 



P(B fw > x(l - p)), x^oo, (2. 



where B- m is the service time of the m-th customer (counting backwards) in the elapsed busy 
period. 

Proof. Let us concentrate on the residual busy period as seen by an arbitrary customer 
( "customer 0" ) arriving at time To = 0. With V- m we denote the amount of work in the system 
found by customer — m and by Z_ m the consecutive busy period if V- m = 0. Furthermore, 
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for m > 0, T_ m is the time between the arrival of customer —m and time and T m is the 
time of arrival of the m-th customer after time 0. We may write 

oo 

P(Z rp >x) = ^P(y_ m = 0,Z_ m >T_ m + zr) 

m=l 
oo 

= P ( V ~™ = °) P ( Z -m > T -m + X) 

m=l 
^ x 

= =- V P(Z > T m + X). 

br *■ — ' 

Prom this, the proof is quite straightforward in the case of constant inter-arrival times T m = 
ma. In that case it follows from Proposition 2.1, that for any 5 £ (0, 1) and d± > 1 there is 
an xq such that 

P{Z > T m + x) = P(Z >ma + x)< (1 + 5)EtP (B > di(x + ma)(l - p)) , 
for all x > xq and m > 1. For x > xo this gives 

oo 

P(Z rp >x) < (1 + 5) P(B>d 1 (x + ma)(l-p)) 

m=l 



Now let 5 — > 0, c?i — > 1 and use that £?^™ € Z7£V (by Property (7e) in Appendix [A[) to obtain 
the desired upper bound 



P{Z rp >x) < (l + o(l))^-P(V w >x(l-p)^j , x^ 



oo. 



The lower bound can be derived similarly. 

When inter-arrival times are not constant the proof is more involved since Zq and T m are 
not independent. First we note that since B € C, then, for any e > 0, 

e - ex = o (P (B > x)) , x^oo. (2.9) 

We shall now develop upper and lower bounds for ^m=o P (^o > T m + x), which coincide for 
x — > oo. 

Upper Bound. For any e £ (0, 1), 

P(Z > T m + x) < P (Z Q > -ex + ma(l - e) + x) + P(T m < -ex + ma(l - e)). 



From Proposition 2.1 for any d\ £ (0, 1), 

oo oo 

Y P(Z > x{l - e) + ma(l - e)) < (1 + o(l))Er ^ P(B > d x (l - e)(l - p)(x + ma)), 



m=0 m=0 



as x —7- oo. For notational convenience we set c\ = d\(l — e) and note that c\ f 1 when di f 1 
and e \. 0. Furthermore, 

oo 

VP(B> Cl (l-p)(x + ma)) ~ -— P—-P(B* W >c x {l-p)x). 

^(1-p) 



5 



We use t n to denote the inter-arrival time of customer n and customer n + 1, thus, T m = 
t\ H + t m . By the Chernoff inequality we have, for any r > 0, 

P(T m < -ex + m«(l - e)) = P (V rTm > e "s=-rma(i- e )^ < (E e - ril ) m e ma(1 - E) - rc3; . 

Since Eti = a and e > 0, we can choose r > sufficiently small, such that 

e ™(i-s) Ee -rti < L (2.10) 

Then 



E P ( T ™ ^ ~ £X + m «( X - £ )) ^ e ~ rCT £ (e ra(1 " e) Ee- rtl ) m = - 



m=0 



m=0 



e ra(l-e)^ e -rii 



• (2-11) 



Thus, from (2.9), 



Em=0 P ( Z > T m + 



X) 



lim sup ^ r , - ■_, 

*->5o P(S/« >ci(l-p)x) " ci(l-p) 



< Er- 



Since £?-^ w G Z1ZV (Property (7e) in Appendix |A[), letting c\ to 1, we get 

,■ Em=o P ( Z o >T m + x) p 
lim sup — £7 „ j... - ; — ; — < Er 



P(Bf w > (1 - p)x) 



1-p' 



which concludes the upper bound. 

Lower Bound. For any e 6 (0, 1), put 



x(l + e) 



+ em 







where [y\ denotes the integer part of a positive real number y. Obviously, 

P(Z > T m + x) > P(Z > T m+nxm ,T m+nxm -T m >x) 

> P(Z > T m+nx m ) - P(T m+Tlx m -T m <x) 
= P(r > m + n X!m ) - P(T n < x). 



From Proposition 2.1, for any C2 > 1, 

P(r > m + jij,. m ) > (1 + o(l))ErP(B > a(l - p)(m + n Xim )c 2 ). 



Similar to (2.10), we can choose r > such that 

e ra(l+|e)/(l+e)E e -rti < X 

so that 

00 00 00 

-rT„ 



< x) 



P{e~ rT ^ > e" ra ) < e rx Ee~ rT ^™ = ^ e rx (Ee~ rtl ) nx 



m=0 



m=0 m=0 
,rx (\z a -rt 1 \ x ( 1 + £ )/ a - 1 



m=0 



< 



< 



e -\erx 



1 - (Ee- r *i) £ - (1 - (Ee- r *i) £ ) Ee~ rt i ' 
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Thus, by (2.9), 



OO p 

2 P(^o > T m + x) > (1 + °( 1 )) C2(1+ P £ w 1 _ p) P (- B/M > C 1 - PM 1 + ^) 



m=0 

Letting C2 ! 1 and e i 0, we get the desired result. 



□ 



Remark 2.1. For the class of RV tails the equivalence (2.6) was proved by Zwart [29] . The 
asymptotic result (2.6) also holds in a class of so-called square-root insensitive subexponential 
distributions under the additional condition that the second moment of the inter-arrival time 
distribution is finite |18j . More precisely, Jelenkovic et al. |18j established the following result 
for the stable Gl/Gl/1 queue. If the following three conditions are satisfied: 

(a) The distribution of service times is square-root insensitive: 

P(B > x + \fx) ~ P(B > x), x -> oo; 

(b) also, the distribution of B belongs to the class of so-called strong concave (SC) distri- 
butions - which is a sub-class of S*; 

(c) the distribution of inter-arrival times has a finite second moment; 



then (2.6) holds. Under the same conditions, it may be shown that the distribution tail of 
the number of customers served in a busy period, r, has similar asymptotics: 

P(t > n) ~ EtP(B > na(l - p)). 

Therefore, one can conclude that the asymptotics (2.7)-(2.8) are also valid under conditions 
(a)-(c) above. It would be worthwhile to formulate and prove Corollaries 2.3 and 2.5 and 
Theorems |2.4| and |3.2| (the main theorem) for the class of service time distributions that 
satisfy (a) and (b) above, under the restriction that the arrival times satisfy (c). 



Proposition 2.2 states that, for large x, the events \J^ = i{B- m > (x + ma)(l — p)} and 
{Z rp > x} are equally likely. In the sequel we shall need that these two events actually occur 
simultaneously (for large x). This statement is made precise in the following corollary. 

Corollary 2.3. IfBe£f]V, then 

oo 

P(Z rp >x) ~ ^P(Z rp >x,5_ m >(x + ma)(l-p)) (2.12) 

m=l 

oo 

~ P(Z rp > x, |J {B- m > (x + ma)(l - p)}), x -> oo. (2.13) 



m=l 
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Proof. First we show that ( |2.12| ) implies ( |2.13[ ). Note that 

oo 

P(Z rp > x, (J {£_ m > (x + ma)(l - p)}) 

m=l 

oo 

< ^ P(Z r P > x, 5_ m > (x + ma)(l - p)) 

m=l 

oo 

< P(Z rp >x,\J {B_ m >(x + ma)(l - p)}) 



m=l 



+ ^ P(-B_ mi > (x + mia)(l - p),B-m 2 > (x + m 2 a)(l - p)), 
where the last sum is not bigger than 



oo oo 



Yl P ( B ~™i >( x + ™ia)(l " p))P(£_ m2 >{x + m 2 a)(l - p)) 

mi=l TO2=1 

~ (pP(S /u, >x)) 2 = o(P(S /w >x)). 
Using Proposition 2.2 this proves that ( |2.12[ ) implies (2.13). 



It remains to show that the right-hand side of (2.12) matches (2.8). As before, we use t n 
to denote the inter-arrival time of customer n and customer n + 1. Assume that for some 
constants e > and R > the following events occur 

1. E%$(x) := <j£- m > (x + R) 1 -P+f( £ 1 +^ + ma(l - p) + ema(l + p) + (1 + e)a + 2i?}; 

2. := {for aU n > 1 : £™ =1 B_ m+ * > n/3(l - e) - i?}; 

3 - ^m? : = { for all n > 1 : Y%=l t-m+i < + e ) + R Y> 
4. := {for all n > 1 : £?=i ^ > na(l - e) - #}; 

then - the amount of work seen upon arrival by customer n - satisfies, for n > — m, 

V n >Yl ~ U+i) >(x + R) 1 ~ P + p) - (n - l)a(l - p + e(l + p)). 
Therefore all customers ra with 

are in the same busy period, so that 



n — 1 < 



x + i? 
a(l — e) ' 



x. 



i=l 



We thus have 



^P(^(x)n^ 2 n^ni?f) < £ P(Z* > s,B > (x + ma)(l -p)) 

m=l 

oo 

< ^P(B>(i + m«)(l-p)), 



m=l 



m=l 
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and we are done if a lower bound for Ylm=i P {^m\{ x ) ^ B m2 ^ \ ^ tfij is shown to 
be arbitrarily close to (1 + o(l))j^-P {B' w > x(l — p)), as x — > oo. For any fixed 5 > 
we can find (by the strong law of large numbers) e and R such that P ( ] > 1 — <5 and 



(E £ ' 



D E% R ) > 1 — S. Thus, as x — > oo, 



E p n ^ n ^ n s 4 *) = E p «5(*)) p «* 2 ) p «* 3 n 5 

m=l m=l 

oo 

>(l-5) 2 EP(<iM 

m=l 

- -P (V >{x + R) l - p ±fg + ^ + (1 + e)o + 2fl) 



l-p + e(l + p) 



l-p + e(l + p) V !- £ 



where we have used B^ w G £. Now, first letting e — > 0, using S^"' G X1ZV (Property (7e) in 



Appendix VM , and then 8 — > the proof is completed. □ 



Remark 2.2. Expression (2.12) shows that the occurrence of a large residual busy period is 
due to a single large service time in the past. This can be explained as follows. The busy 
period is the sum of services of the customers in that busy period. The number of customers in 
the busy period after time (the point of arrival) is almost surely finite. There are, however, 
infinitely many service times in the past, each of them being potentially large. This leads to 
the integrated tail of the service time distribution. 

Besides the busy period and the residual busy period, there is a third entity whose distri- 
bution is the same for all non-idling and non-preemptive service disciplines: B rp , the residual 
service requirement of the customer in service (if any) upon arrival of a new customer. The 
tail asymptotics for the distribution of B rp are determined in the following theorem. Not only 
is the theorem of interest in itself, but several steps in its proof will also be useful in proving 



our main result in Theorem 13.21 

Theorem 2.4. If B £ £f)T> then, for any non-preemptive and non-idling service discipline, 

P {B rp > x) ~ pP (B fw > 



as x —> oo. 



Proof See Appendix [Bj □ 

Corollary 2.5. If B G Cf]T>, then for any non-preemptive and non-idling service discipline, 
the waiting time W and residual busy period Z rp seen by a customer arriving to a stationary 
GI/G/1 queue satisfy W < Z rp a.s. and therefore, as x — >• oo, 

oo 

P(W > x) ~ E p ( w > x i B -m > (x + ma)(l - p)) , x^oo. (2.14) 

m=l 
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Proof. Since the service discipline is non-preemptive we have W > B rp almost surely, so that 



by Theorem 2.4 



P (W > x) > P (B rp > x) ~ pP (B fw > 



x — > oo. 



Thus, in the following, we may neglect terms which are o (P (^B^ w > x)). Using W < Z rp 
(almost surely) and Corollary 2.3 we therefore have (similar to the proof of Theorem 2.4) 

P(W>x) = P (W > x, Z rp > x) 

oo 

~ P(W > x, Z rp > x, (J {B_ m >{x + ma){\ - p)}) 

m=l 

oo 

~ ^P{W > x,Z rp > x,B^ m > {x + ma){l- p)) 

m=l 

oo 

= P ( W > x > B -rn > (x + ma)(l-p)). 



m=l 



□ 



3 Random Order of Service 

We now turn to the GI/GI/1 queue with Random Order of Service. We start with analyzing 
the waiting time conditional on the initial queue length q, none of these customers having 
received previous service. It is convenient to associate service times with customers in their 
order of service instead of their order of arrival. The customer which is served first has a 
service time B\, the second has B2, etc. Denote with Wros(q) the conditional waiting time 
of an arbitrary customer in the queue. 

The following lemma does not require any assumptions on the distributions of service times 
and inter-arrival times. It shows that Wros(q)/q converges in distribution to a random vari- 
able whose distribution has support [0, j^r]. Note that this contrasts with the FCFS queue, 
in which case the corresponding quantity W{q)/q (for the last customer in line) converges to 
the constant (3 almost surely. 

Lemma 3.1. As q —> 00, 

P (w RO s(q) > ^-y\ ((1 - y) + ) ^ , (3.1) 

uniformly in y G [0, 00). 

Proof. Since the limiting distribution is continuous and non-defective, by the monotonicity 
of probability distribution functions, it is sufficient to prove point-wise convergence. 

Let us thus fix y E (0, 1). For n = 1, 2, . . ., denote by Q n the number of customers in the 
queue at the time instant of the nth service completion. We define the event A\ by 

Ai := |q< G [g(l - e) - i(l - p + e), q(l + e) - i(l - p - e)\ V i = 1, 2, . . . , ^— j . 
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By the strong law of large numbers, for any e > 0, there exists q = q{e) such that P(-Ai) > 1— e 
for all q > q. 

Let q~ = q(l — e) — i(l — p + e) and qf = q(l + e) — i(l — p — e) and denote N(v) = 
min{n : Y2i &i > v }l customer N(v) is in service at time v. Defining 



An 



N{v) £ 



-p(l-e)-R, V -{l + e)+R 



0, - y 

. 1 ~P . 



for any e > 0, we may choose R > such that P(^) > 1 — £. Thus, P(Ai f"| A2) > 1 — 2e and 



P Wkoste) > 



i-p< 



where 



We further define u 

fa 



P{y) := P (V ROS (<z) > Y^- p V:M[\A^j 



|(1 - e) - R and n+(«) = |(1 + e) + R. Since 



{Wros(<7) > j^y} implies that customer was not selected in the first N ( izr^y) trials, we 
have, as q and u tend to infinity keeping y constant, 



P{y) < n^ (u) ( 1 



1 



I n (u) 
(l + o(l))exp - ~ 



(l + (l))exp^-^ 
(l + o(l)) (l-i/ 



n (it) 



1 



g(l + e) -u(l-p-e 
(1 - e)(l - p- £)\h^ 



(l + e)(l-p) 
(1 + o(l))(l - y + O(e))^ . 



Similarly, 



P(y) > Utt ] ( 1 



1 



0(e) = (1 + o(l)) (1 - y - 0{e))*=p=i - 0{e). 



Letting e pass to 0, we obtain (3.1). 



□ 



The main result of our paper is stated in the next theorem. 
Theorem 3.2. In the GI/G/1 ROS queue with B G Cf]T>, we have 
P(W RO s>x) ~ pP(B^>x) 

fva+x 

dv I dP(B < z 



+ 



+ 



+ 



dv 



CI 

00 



(va+x)(l— p) 
va+x 

dP{B <z)[l 



CO 

va 



j A _ (x + ^a-z)(l-p) ^ ^ 

1 

{x + va — z)(l — p) \ !-p 



dP(5 < z ) 1 - 



pz 

x(l-p) 



1 

l-p 
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where c 



ap ' 

Proof. In Appendix [C] □ 



Letting W^ os be a random variable independent of B with the limiting distribution of 
Wnos(q) / q, as q — > oo, we may conveniently rewrite the above as: 

P(W ROS >x) ~ P(^>x) 

+ dvP ((va + x)(l-p)<B<va + x, W& os > a ( x + v ® ~ 

+ J°° dvp(va<B<va + x,W£ os > a ( x + v ®- B ) ^ 

+ £ A,P + *)(!- p) < B < «a, WSos > g _i(!_ p) ) • 

This allows for the following interpretation: The waiting time is larger than x when one of 
the following occurs: 

1. (first term) The customer in service has a residual service time larger than x. Recall 
that, by Theorem^ P (B r P > x) ~ pP (B? w > x). 



2. (second term) A customer (with index —v) that arrived at some time —t = —av between 
time — acx and time required a service z larger than (t+x)(l—p) but smaller than t+x. 
The service times of other customers in the system at time —t are negligible compared 
to z. The large service time ends at time — t + z £ (0, x), leaving approximately z/a 
competing customers in the system. Customer thus waits for approximately — t + z + 
W^ os z/a. Thus Wp^os needs to be larger than (x + t — z)a/z. 

3. (third term) This term is similar to the previous. Now, the large customer arrived at 
time —t < —acx with a service requirement z S (t,t + x). That customer thus leaves 
at time —t + z£ (0, x) with approximately z/a customers in the system. 

4. (fourth term) Again, the large customer arrived at time —t < —acx, but leaves before 
time 0: its service requirement is z £ ((t + x)(l — p), t). Neglecting the size of the 
customer in service at time 0, the "service lottery" starts immediately upon arrival of 
customer 0. The number of competing customers at time is approximately t/a — (t — 
z)//3 which is the number of arrivals minus the number of departures between times —t 
and 0. Therefore, the waiting time of customer is larger than x if VFj^Qg is larger 
than x/(t/a - (t - z)/(3) = f3x/(z - t(l - p)). 

4 The M/G/l queue with regularly varying service time dis- 
tribution 

In this section we restrict ourselves to the case of Poisson arrivals and regularly varying service 
time distribution. In this case we are able to obtain detailed tail asymptotics for the waiting 



time distribution by applying Lemma A.l for Laplace-Stieltjes transforms (LST) of regularly 
varying distributions to an expression of Le Gall [22] for the waiting time LST in the M/G/l 
queue with ROS. In Section [5j we shall present an alternative approach to the same result, 
viz., we shall work out the general tail asymptotics of Section [3] for this case. 
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We consider an M/G/l queue with arrival rate A = 1/a and service time distribution 
B(-) with mean (3 and LST As before, the load of the queue is p = A/3 < 1. 

The LST of the waiting time distribution for the ROS discipline is given by (see Le Gall [22 
or Cohen [T3], p. 439): 

E[ e -WHos] = l- p + p ^{ s ) - f ^( S , Z )V( S , z)z, (4.1) 



with 



* (1 _ z) ^ + A(l-z)}-/3{A(l-z)} 



exp 



z-f3{X(l-z)} 

y 



y _ p {s + A (l - y)} 
where (3^ w {-} is the LST of the forward recurrence time of the service time: 



(4.2) 
(4.3) 



P fw {s} 



def 1 - P{s} 



and p{s] is the LST of the busy period distribution, satisfying the relation 

p{s}-(3{s + X(l-p{s})} =0. 



(4.4) 



(4.5) 



It is possible to rewrite (4.1) in a simpler form using integration by parts. Indeed 



/j,{s} 



dz 



(s,z)^(s,z)z 



®(s,z)^(s,z) 



1 



<£(s, z)^(s, z)z 



and the simple relations 



®{s,p(s)) = 1 - ti{s}, 



J^{s} Z- P{s + \(l- Z)}' 

1 - (3{s} 



l-p ' 



yield 



*(s,p{s}) = 0, *(s,l) = l, 
1-P{s} 



&(s,z)i$>(s,z) 



p{s} 



l-p 



Using this relation and ( |4.4[ ) in (4.1) yields the following simpler alternative: 

P(l-P) f 1 



with 



$(s,z) 



1 + 



<$>(s,z) 



&(s, z)^(s, z)z, 



(4.6) 



l-p 



z - p{s + A(l - z)} 



1-z 



_§s_ 
l-p 



1 - Z 



z - p{s + A(l - z)} z-p{\(l-z)}' 
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As before, we write f(x) ~ g(x) if f(x)/g(x) — > 1 when x —> oo, and similarly we write 
a(s) ~ b(s) if a(s)/b(s) — > 1 when s — )■ 0. In this section and the next we assume that the 
service requirement distribution B(-) is regularly varying with index — v, 1 < v < 2: 

C 

P(B > x) ~ — — -rx~ v L{x), x ^ oo, (4.7) 
with C a constant, T(-) the Gamma function and L(-) a slowly varying function at infinity, 

cf. |5]. 

Lemma A.l in Appendix [A| implies, in combination with the assumption that (4.7) holds 
with v E (1,2): 

/3{s} - 1 + /3s ~ CVL(l/s), as s -> 0. (4.8) 



In addition, cf. (4.4), 

1 - P fw {s\ ~ ^V-^fl/a), as s -> 0. (4.9) 

/? 

De Meyer and Teugels [23j have proven that the busy period distribution of an M/G/l queue 

with regularly varying service time distribution is also regularly varying at infinity. More 
precisely: 

MM - 1 + ~ _ C s v L(l/s), as s -> 0. (4.10) 

We shall use the first-order behavior of /x{s} further on. 
Our goal is to prove the following theorem. 

Theorem 4.1. // the service time distribution in the M/G/l queue operating under the ROS 
discipline is regularly varying at infinity with index —v £ (—2,-1), then the waiting time 
distribution is regularly varying at infinity with index 1 — v £ ( — 1, 0). More precisely, if (4-7) 
holds then, as x —> oo, 

P(Wkos>x) ~ T £ —h{p,v)P{Bf v >>x) 
l-p 



l-p ^ 'T{2-v)P 



with 



h{p,u) = [ f(u,p,v)u, (4.11) 

J 

Hu.r.v) * T ^( i ^)"" 1 (1 -„,ii ; + ( 1 + I ^ ) )''(i-„)^. (4.12) 

Remark 4.1. The following result for the waiting time distribution in the M/G/l queue 
operating under the FCFS discipline is well-known (cf. Cohen [12] for the regularly varying 
case; see Pakes [23] for an extension to the larger class of subexponential residual service time 
distributions): if (4.7) holds, then 

P(Wfcfs>x) —P(B fw >x), x^oo. (4.13) 

l-p 

We can now conclude that 

P(Wrqs > x) ~ h(p, i/)P(Wfcfs > a), a; -»• oo. (4.14) 
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Figure 1: A plot of the function (p, v) h-> h(v, p) for < p < 1 and 1 < v < 2. Note that the 
minimal value for h seems to be h(l, 3/2) = T(3/2) ~ 0.88622 



Remark 4.2. It is actually possible to prove that h(p,is) is always less than 1: indeed, the 
function f(u, p, v) is strictly convex in v (as a sum of exponentials) and therefore h(p,u) is 
also strictly convex in v\ moreover, simple computations using integration by parts show that 



/ f(u,pA)V-= f(u,p,2)u = l, 
Jo Jo 



and therefore 

h(p, v) < 1, for all 1< v < 2, < p < 1. (4.15) 

Numerical computations with MAPLE suggest that h(p, u) is decreasing in p and thus 
always larger than its limit when p — > 1, which is by simple arguments equal to Un- 
fortunately, we have not been able to find a simple proof for this fact. In Figure 1 we have 
plotted h(p, v) for < p < 1 and 1 < v < 2. 

It is interesting to observe that the tail behavior of Wros and Wfcfs is so similar in the regu- 
larly varying case. This strongly contrasts with the situation for the M/M/l queue, where the 
purely exponential waiting time tail for FCFS strongly deviates from the Cox~ 5 ^ e e~ ClX ~ C2xl 
waiting time tail behavior that was exposed by Flatto pTSJ for ROS. 

Remark 4.3. The first part of h(p, u) is a Beta function, and the second part is a hypergeo- 
metric function (cf. Abramowitz and Stegun [1]). In particular, 



Jo 



i / i \ r(i/)r(^- + i) 



i-p > r(y + I ^ + i 



Using partial integration and the above formula for Beta functions, one gets the following 
form, which is useful for future comparisons (see Section [5]): 

KPiv) = l-p+l g(u,p,v)du, (4.16) 



with 



g{u,p,u) = 



i-pj V 1 ~p 



i 

ii-p. 
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Proof of Theorem 4-1 ■ We shall prove the theorem by applying Lemma |A.1| to an expression 
for the LST of the waiting time distribution. So we need to consider E[e _sH/ROS ] — 1 as s — > 0. 
In order to do that, we use the change of variable 1 — z = e(s)u, with e(s) = 1 — p{s}, in 
(4.6) to obtain 



fj,{s} 



®(s, z)^f(s, z)z = £ 



(s) 

Jo 



s, 1 — e(s)u)^>(s, 1 — e(s)u)u. 



Let us first evaluate the function ^. One can write 



log*(s,l-e(s)«) = - / - 

Jo 1 



e(s)uv 



e(s)uv — j3{s + Xe(s)uv} ' 



and, using (4.8) and (4.10), 



lim 



1 



1 — e(s)uv — (3{s + Xe(s)uv} = (1 — p){\ — uv), 



s->o e(s) 

the limit being uniform in u and v. Therefore the following limit holds uniformly in u: 



lim ^(s, 1 — e(s)u) = exp 
s— >o 



1 



uv 



1 — p Jq 1 — UV 



(1 - u)i-p. 



(4.17) 



The evaluation of $ is not difficult either. First note that the denominators appearing in 
1 — e(s)u) can be expressed in terms of @* w {-} as follows: 

1 — e(s)u — (3{s + Xe(s)u} = 

-(1 - p)e{s)u + (3s- (/3s + pe{s)u)(l - (3 fw {s + Xs(s)u}), 
1 — e(s)u — f3{Xe(s)u} = 

-(1 - p)e{s)u - pe(s)u(l - (3 fw {Xe(s)u}). 



Then write 



(3s \ 

:(s)u ) (1 — e(s)u — (3{Xe(s)u}) — e(s)u(l — e(s)u — (3{s + Xe(s)u}) 

1-pJ 

= e{s)u\((3s + pe{s)u)(l - f3 fw {s + Xe(s)u}) 



+ 



_pJ3s_ 
1-p 



pe(s)u){l - p fw {Xe(s)u}) 



and finally, as s — > 0, 

$(s, 1 - e(s)u) 



1 



(l-p)(l-u) 



+ 



pu 



p(l-u) 



(1 - f3f w {s + Xe(s)u}) 



1-p 



(1 - (3^{Xe(s)u}) 



that 



We take into account now the regular variation assumption (4.7), which yields (4.9) so 

Cs v - X 



$(s, 1 — e(s)u) ~ 



p(l - u) ( pu / 1 



(3(l-p)(l-u) 1 1-p VI -p 

pit 



Ae(s)« 



+ 1 + 



1-p 



1 

s + Ae(s)w 
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Using Potter's Theorem (see (23], Theorem 1.5.6), setting 5 = {y — l)/2, there exists X > 
such that, as long as s < 1/X and Xe(s) < 1/X, 



L 



1 



L 



Xe(s)u 
1 



< 2 max 

< 2 max 



-6 



Xe(s)uJ \Xe(s)u 
1 \ s ( 1 



1 + ^u 



L(l/s) 
6 



l + ^Mu 



L(l/s). 



+ Ae(s)u/ 

These bounds allow application of the Dominated Convergence Theorem to 

r 1 i 

which yields 

1 - E[c- sH/ROS ] ~ / /(«, p, i/)u s^ 1 L(l/s), as s 0. 

I — P L./n 



(4.18) 
□ 



Using Lemma |A.1[ the theorem follows. 

5 Agreement of results 

While Sections [3] and [4] use completely different methods of proof, it is clear that the asymp- 
totics for the M/G/l queue with regularly varying service time distribution (as in Theo- 



rem 



4.1 ) has to be a mere consequence of Theorem 3.2 This section shows that this is indeed 



true. As a first step, we give another asymptotic expression for P(W / ros > x ) which, while less 
intuitive than the rewriting proposed in Section [3j bears a strong similarity with Theorem 4.1 



Lemma 5.1. In the GI/G/1 ROS queue with B G C{~}T), we have 



P(^ros>x) ~ P P(Bf w >x) 



+ 



1 



o \a(l-p) 

x 



Bl x(1 P ) <b< *{i-p) } 

L mi.-\- \ — n — mi, 1 



+ 



P B > 



pu+l — p 

xQ- - p) 
pu 



1 — u) du 



Proof. To simplify the computations, assume that B has a density function b. This is, how- 
ever, not really needed for the result. 



The first term in Theorem 3.2 coincides with that in the expression above. Next, consider 



the second and third terms in Theorem 3.2 together, using the changes of variable z \— > u 
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(x + v a — z)(l — p)/(pz) and v ^ t = (x + va)(l — p)/(pu + 1 — p). 

!" *, r + ' d , b(z) (i - * 

J J(va+x)(l-p) V P z J 

- r*r d»K, ) ('i-t' + '"-^ i -^ A 

io y(«a+x)fl-p) V P z ) 





ca; J (ua+a:)(l — p) V P^ 



oo />! 



(x + f a)(l — p)\ p(x + va)(l — p) 



' 1 pu + 1 - p J (pu + 1 - p) 2 K ' 



J 

~d„ / f d,6 (I* + g^Hl -P) (1 _ M) A 

7 v pu+i-p ) {pu + i- P y 

1 du (1 - u)^ P dt> 6 ^ + " ^ P(X + " a)(1 " P) 

JO 

1 



d«(l-«)i-p / dtt&(t). 

pu+1— p 



pu + 1 — p J (pu + 1- p) 2 



Finally, focus on the last term of Theorem |3.2| We use the changes of variables z ^ w 
z — va(l — p) and then w ^ u = x(l — p)/w. 



oo pva 



dv I dzb(z)(l X(1 P) ^ " 



cx J(va+x)(l-p) \ Z-Va(l-p) 

r dv r dw b(w + vail - p)\ ( 1 - X(1 ~ A V"" 

■/<» A'(i-p) V u> / 

/°° dw fl - X(1 1 " P f°° dvb(w + va(l- p)) 

1 Z" 00 / z(l-p)\T^P „ / 10 

' dw 1 - — P[B> 



«(1 - p) Jx{i- P ) V w 7 V P. 

= X - I" \?(B> X ^-\{l-u)^du. (5.1) 
a Jo u 2 \ pu J 

The proof of the lemma is completed by collecting the terms. □ 
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Remark 5.1. It is interesting to note that (5.1) can be rewritten as follows: 
-1 



.r 



1 ;P(b> X{1 P) ) (1-U)^du 



a ./q u \ P u 

x 
a 

x p 



du 



, P(B>z)p(w^ os > — )dz 
a x(l - p) J *a-p) V KUb zJ 

p 

p 2 _ p ( Bfw > ECW) iB / ttwSog > ax ) m 



1 - p \ P 

In the case of Poisson arrivals and regularly varying service times, assuming again that 

C 



P (B > x) 



r(i-iz) 



x~ v L(x), 



and recalling that A = 1/a, it is easy to go from the expression in Lemma 5.1 to the expression 

CL(x)x l - u 



P(W ROS > x) 



a{l-p)Y(2-p) 



l 

1 - P + I a(u,p,^)du 





Using (4.16) it is seen that this corresponds to Theorem 4.1 



6 Heavy-traffic limit for the waiting time distribution 

In this section, we consider the M/G/l queue with general service time distribution. We are 
interested in the heavy-traffic case: p — > 1. The main result is to find a sequence A(p) which 
tends to as p — > 1 such that E[e-" A ^) M/ROS ] tends to a proper limit as p — > 1. This way we 
shall be able to retrieve a result of Kingman [20] for the case of finite service time variance, 
as well as derive a new result for the case of regularly varying service time distribution with 
infinite variance. 

The following lemma will be useful in the sequel. 



Lemma 6.1. The following bound holds for all s > and < z < 1: 

1- an 



^(s, z) < exp 



/3s 



Moreover, for any t > ; 



lim 1 - fist) = e~*. 

(p,«H(i.o) 



Proof. Using the inequality fi{s} > 1 — /3s, one finds 

y - (3{s + A(l - y)} < y - 1 + /3s + p(l - y) < /3s, 



and (6.1) follows from the definition (4.3), since 



log^(s,z) > 



1 dy _ 1 - z 
Js ~ /3s ' 



(6.1) 
(6.2) 
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The proof of (6.2) follows a similar argument: indeed, 

y 



log*(s, 1 - Pat) = 



fist 



Pstv 



y - P{s + A(l - y)} J 1- Pstv- P{s + pstv}' 



and, as s — > 0, 



1 - Pstv - P{s + pstv} ~ Ps(l - (1 - p)tv). 



This yields (6.2) and concludes the proof of the lemma. 



□ 



The LST of the steady-state waiting time distribution under FCFS is given by (cf. [13] , p. 
255): for any s > 0, 



£J e -sW FC FS] 



1-P 



1 - ppf w {s}' 

The following lemma illustrates the relation between Wros an d Wfcfs * n heavy traffic. 

Lemma 6.2. Assume that P < oo and that there exists A(p) > which can serve as a proper 
scaling for Wfcfs-' for any oj > 0, 



lim E[e _ajA( - p ^ FCFS ] = E[e~ wM/FCFS ], 



(6.3) 



where Wfcfs is a non-negative random variable. Then A(p) is a proper scaling for Wros-' 
for any to > 0, 

poo - - 

lim E[e _wA< - p ^ w/ROS ] = / E[e~ w ' M/FCFS ]e~*i = f E[e~ wW/ROS ]. 
p^ 1 Jo 



Proof. The starting point of the proof is Equation (|4.6|). First, using integration by parts, 

1-z 



1±_ 

1-P 



M { S } z - p{s + A(l - z)} 



W(s,z)z 



Ps 



+ 



1 - P J^s} 



*(s,«)z, 



and the integral above can be bounded using (6.1) as 

■i 



H{s} 



\I/(s, z)z < I exp 

H{s} 



1 - Z 



z < /3s. 



The second part of (4.6) can be expressed in terms of Wfcfs : 

-^(s,z)z 



{s} z - P{X(1 - z)} 

pi | 

= - / E[c- a ( 1 -^ h/fcfs ]^(s,z)z 

JaU} 1 - P 



[ 88 E[e- pstWFCFS ]^(s,l- Pst)t. 
1 — P Jo 



Plugging these two relations into (4.6) yields 
E [ e -sty RO s] = p 



£(£) 
/3s 



E[e-P stWFCFS }^(s, 1 - pst)t + 0{l - p), 



(6.4) 
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uniformly in s > 0. We now show that 



ft eJ ^r - *»■ (6 ' 5 » 

with, as before, e(s) = 1 — p{s} and /i{s} is the LST of the busy period distribution. Let 
fip(s) be the LST of the busy period of the M/G/l with service time distributions £?(•) and 
load p. Obviously, for p > p we have p(s) < pp(s), for all s > 0. Hence, 

liminf ^^ >lim l-^A(,))_ fa 



P-+1 A(p) ~p->i A(p) 1 



This is true for any fixed /) 6 (0, 1). Letting /? pass to 1 we obtain (6.5). 



Finally, under Assumption (6.3), Lemma 6.1 allows to apply the Dominated Convergence 



Theorem to (6.4), and the lemma is proved. □ 



Using Feller's continuity theorem, Lemma |6.2| can be rewritten in a more compelling way. 

Corollary 6.3. Assume that there exists A(p) > 0, and a random variable Wfcfs> such that 
the following limit holds in distribution: 

lim A(p)Wfcfs = Wfcfs- 



Then, in distribution, 



lim A(p)Wros = ^WWs = Wros, (6.6) 



where Y is an exponential random variable with mean 1, independent ofWpcFS- 

Remark 6.1. In vi ew o f the PASTA property and the fact that the workload is the same under 



FCFS and ROS, (6.6) states that the scaled waiting time Wros equals in distribution the 
product of the unit exponential Y and the scaled workload Vros = Wfcfs • Put differently, 
Pl^ROS > x|Vros = y) = P(X > x /y) = e~ x / y . The latter result might be intuitively 
understood by referring to a snapshot principle. Consider a tagged customer. If it is not 
being elected for service at the end of some ROS services, then in the heavy-traffic limit the 
remaining workload that it sees has not changed. The randomness ('memory less') property 
of ROS then implies that the remaining waiting time of the tagged customer has the same 
distribution as before. 



As a first application of Lemma 6.2, consider the case where the service time distribution 
has a finite second moment. The following theorem extends a result of Kingman [2D], where 
it was additionally assumed that (3{s} exists for some s < 0. 

Theorem 6.4. Assume that the variance a 2 of the service time is finite and let 



Then, for any u > 0, 

lim E[e- a;A ^ iyROS l = / 
p^ 1 Jo 



oo e -t 



1+ojV 
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Proof. It is known from a classical result of Kingman [19] that, with A(p) defined by (6.7), 

lim E[e- wA ^ iyFCFS l = , 



so that we may apply Lemma 6.2 □ 



Remark 6.2. Thus, when the service times have a finite variance, Wfcfs has an exponential 
distribution, so that Wros is the product of two independent exponentials with unit mean. 
Letting K\(-) be the modified Bessel function of the second kind, simple calculations show 
that 



/>oo 

P(WrOS > X) = / P(iW FC FS > x)e"*t 
J 



X 

1 

t 



t 



exp 

which coincides with Theorem 6 of |20j . 

In the case of a service time distribution with regularly varying tail and infinite variance 
a similar, but new, result can be obtained from results of Boxma and Cohen [8j for FCFS. 



Theorem 6.5. Under Assumption (4-8), with 1 < u < 2, let A(p) be the unique root of the 
equation 

\Cx v - x L{x) = 1 - p, x>0, (6.8) 
such that A(p) 1 for p f 1. Then 

poo —t 
lim E[e- wA(p)WROS ] = / j . T t 

Proof. It has been proved in [8] that A( / o) exists and that for any oj > 0, under Assump- 
tion ( gig) ), 

lim E[e- wA( ^^ FCPS ] = -. 

The theorem now follows from Lemma 16^21 □ 



Appendices 

A Classes of distributions 

Definitions and properties 

We say that a random variable belongs to a certain class if its distribution function belongs 
to that class. 

1. A cdf F belongs to the class C of long-tailed distributions if there exists a y > (or, 
equivalently, for all y > 0) such that, as x — > oo, 

F(x + y) ) 1 
F(x) 
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(a) If F G £, c > 0, and G is another distribution such that G{x) ~ cF(x) as x -> oo, 
then G G £. 

(b) If F G £ and m + = m + (F) = J °° F(t)dt is finite, then the integrated tail distri- 
bution F 1 belongs to £ too, but the converse is not true, in general. Here 



F 7 (x) = max^0;l-^ F(t)dtj . 



(c) F 1 G £ if and only if F(x) = o(F I (x)) as x — ?■ oo. 



2. A cdf F belongs to the class 1ZV of regularly varying distributions if there exists a v > 
such that 

F(x) = x~ v L{x), 

where L(x) is a slowly varying (at infinity) function. 

3. A cdf F belongs to the class V if 

r F(2X) 

inf ^ — '- > 0. 
x>0 F(x) 

(a) If F G V and m+ < oo, then F J G V. 

4. A cdf F belongs to the class I1ZV of intermediate regularly varying distributions if 

limliminf - = 1. 

5. A cdf F on the positive half- line belongs to the class S of subexponential distributions 
if 



f F(dt)F(x - t) ~ F(x) as x -> oo. 



(a) A cdf F on the real line belongs to S if F(x)I(x > 0) belongs to S. 
6. A cdf F belongs to the class S* if m + is finite and 



rx 

\ F(y)F(x — y)dy ~ 2m + F(x) as x — > oo. 



7. Relations 



(a) [14, p. 50] TZV C ZTZV cC{~}VcS. 

(b) [21] If F G 5*, then F G 5 and F J G 5. 

(c) [2T] If F G £ f] ^ and if m+ is finite > then F G 5 *' 

(d) If F G P and F has an eventually non-increasing density, then F G T1ZV. Indeed, 
put 

K = sup _ . 

x >o F(x) 
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Then, for £ > 1 and for all sufficiently large x, 



F(x) - F((x) _ F{x/2) F(x) - F((x 
F(x) F(x) F(x/2) 

tm dt <K (c-i)xf(x) 



< K 



\xf{x) 



2K(( - 1) ->• 



as C 4 1 • 



(e) It follows from Properties (3a) and (7d) that if F E T> and m + < oo, then F 1 G 
XUV. 



There exists a very useful relation between the tail behavior of a regularly varying probability 
distribution and the behavior of its LST near the origin. That relation often enables one to 
conclude from the form of the LST of a distribution, that the distribution itself is regularly 
varying at infinity. We present this relation in Lemma |A.l below. We use this in Section [4] 
to prove that the waiting time distribution in the M/G/l queue under the ROS discipline is 
regularly varying at infinity if the service time distribution is regularly varying at infinity. 

Let .F(-) be the distribution of a non-negative random variable, with LST <p{s} and finite 
first n moments [ii,...,(J, n (and /iq = 1). Define 



i{s} |S L {s} _^. 



j=0 



Lemma A.l. Let n < v < n+ 1, C > 0. The following statements are equivalent: 

(j) n {s} = (C + o(l))s"L(l/s), 4 0, s real, 



l-F(x) = (C + o(l)) r | ■ ' x-"L{x), x^oo. 



The case C > is due to Bingham and Doney [1]. The case C = was first obtained 
by Vincent Dumas, and is treated in |10j . Lemma 2.2. The case of an integer v is more 
complicated; see Theorem 8.1.6 and Chapter 3 of [S]. 

B Proof of Theorem 12.4 

Note that the distribution of the residual service time of the customer in service is the same 
for all non- preemptive and non-idling service disciplines. We may therefore concentrate on 
the FCFS discipline. 

As before, V- n denotes the amount of work in the system upon arrival of customer — n 
and T_ n is the time between arrival of customer — n and time (which is the arrival time of 
customer 0). In the sequel the random variable V has the stationary workload distribution. 

Lower bound. For any e G (0,1), choose K\ > such that P(V > Ki) < e. Then 
choose an integer n > such that P(T_„ > K\) > 1 — e. Third, choose K<i > K\ such that 
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P(T_„ G (Ki,K 2 )) > 1 - 2e. Then, as x -> oo, 

oo 

P(B rp >x) > P(V- m < Kx,T. n e (K 1 ,Kz),B- m >x + T- m ) 

m=n 

oo 

^ E P ( y -'» ^ ^i. T -n G (K 1 ,K 2 ),B- m > x + T_ m - T_„ + Xa) 

oo 

> (l-2e)(l-e)^P(B >x + K 2 + r_i) 

z=o 

> (l-2e)(l-e)pP(5^>x + iir 2 ) 
~ (l-2e)(l-e) P P(B fw > x). 

Letting e — > 0, we get the correct lower bound. Since B S V, the lower bound is of order 
0{P(Zf w > cx)) for any positive c. 

Upper bound. Let y be a positive number and r](y) = min{n > 1 : Y27=l > y}i 
x(y) = X/i Bi — y. Since is finite, it follows from basic renewal theory that the family of 
distributions of random variables {x(y), V > 0} is tight, i.e. u(x) = sup y>0 P(x(y) > x) — >■ 
as x — > oo. 

Since l? rp < Z rp almost surely, we have by Corollary |2 . 3| and by the lower bound obtained, 
P(B rp >x) = P {B Tp > x, Z rp > x) 

oo 

= (1 + o{l))P(B rp > x, Z rp >x,\J {B_ m >{x + ma)(l - p)}) 

m=l 

oo 

= (1 + o(l)) ^ P(5 rp > x, Z rp > x, £_ m > (x + ma){\ - p)). 

m=l 

Denote by f m (x) the m-th term in the latter sum. For any e G (0, 1), choose K > such that 
P(V > K) < e. Then 

f m (x) < P{V- m <K,B_ m >x + T- m -K) 

+ P(V_ m > K, S_ m > (x + ma)(l — p)) 

+ P(B rp > x,V- m < K, B^ m G ((x + ma)(l - p) - K,x + T^ m — K); 31 < I < m 

+ E / P(V„ m < K, V- m + B_ m G dt)P( X (x + T_ m -t)> X | T— m) 

J (i+ma)(l-p) 

= fm,l{x) + fm,2{x) + f m Ax) + f m Ax). 



Here 






fm,l{ 




< 


/m,2< 




< 


/m,3( 




< 






< 


fm,i{ 




< 



TO— 1 

P(Z rp > X 



p , r _ ()) P(-B_ m > (x + ma)(l-p)); 
fmA x ) ^ u(j;)P(B- m >(i + ma)(l-/))). 
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Since u(x) — > and P(Z rp > x) — > as x — > oo, 

oo 

Y Um,2{x) + f m , 3 (x) + f mA (x)) < (1 + o(l))eP(Z r P > X). 



m=l 

Further, 



Y fm,l{x) < (1 + o(l))pP(fl' tt > x - K) ~ pP(B^ > x). 



m=l 



Letting e 1 0, the proof is completed. □ 



C Proof of Theorem 3.2 



We focus on the waiting time VFros of customer arriving at time 0. Our proof consists of 
three main parts, each corresponding to a typical scenario in which the large waiting time 



arises. The intuition behind these typical scenarios was discussed below Theorem 3.2 (it is 
convenient to treat the two "middle terms" as one scenario). 

Before proceeding, we note that the distribution of the waiting time of customer is not 
affected if we choose to use the FCFS discipline before time and the ROS discipline after 
time 0. Thus, VFros = Wros> wnere Wros denotes the waiting time of customer under the 
modified service discipline. 



The starting point of the proof is (2.14) in Corollary 2.5, which we repeat for convenience 



(the modified service discipline is non-idling and non- preemptive) : 

oo 

p (Wros > x) ~ Y P ( W kos > x,B- m >{x + roa)(l - pj) , x -> oo. 

m=l 

In the verbal discussion, we interpret this relation as follows: there is one "large customer", 
i.e., customer — m for which £?_ m > (x + ma)(l — p), that causes the large waiting time of 
customer 0. Note that any scenario in which the service of this large customer did not start 
before time may be neglected (i.e., is of the order o{P{Bf w > x))): 

oo 

Y P ( W *os > x, >(i + ma)(l- p),T_ m < y_ m ) 

m=l 

oo 

< y p ( B - m >( x+ ma )( x - />)> T - m ^ v -™) 

771=1 
oo 

< Y p ( B -m >(x + ma)(i - P ), r_ m < v- m , y_ m > k) 

771=1 

OO 

+ Y P ( B ~m >(* + " P),T-m < V- m < K) 

777=1 

OO 

< P{V > K) Y P ( B -m >(x + ma)(l - p)) 



777=1 

OO 



+ P(T_ M < K) Y P ( B -m >(% + ma)(l - p)) + o{P{B^ > x)) 

m>M 

(P(V >K) + P(T_ M < K)) -f—P{B^ > x), 

1 - p 
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which may be neglected after first taking M — > oo and then K — Y oo. We have used that 
the sum of a finite number of terms in the above summations is of the order o(P(B^ w > x)), 
see Property (lc) in Appendix [Aj This property, as well as other steps taken in the proof of 
Theorem 2.4 will be used frequently in the following. We have thus proved that, as x — > oo, 

oo 

P(^ros > ~ E P ( W *os > ^ V-m > (x + ma)(l - p), V- < T_ m ) + o(P(B fw > x)). 

m=l 

Part I. We start with the scenario that the large customer is in service for the entire interval 
(0, x). This is the case when the workload V- m < T- m and B- m > T_ m — V- m + x. This 
immediately implies that the waiting time of customer exceeds x. By Theorem 2.4 we have: 

oo 

P (Waos > x, B- m > (x + ma)(l - p),V- m T— m , B— m > T— m V— m ~i~ X) 



m=l 

oo 



£ P (#-m > (x + ma)(l - p), V_ 



m=l 



P(B rp > x). 



This corresponds to the first term in Theorem 3.2 



Part II. We now investigate the event that W^ os > x occurs while the large customer 
(customer — m) is still in service at time 0, but not anymore at time x. Thus, V- m < T_ m < 
B-m + V- m < T- m + x. Let W^ os (q) be the remaining waiting time of customer after 
the first service completion after time if q is the number of competing customers at that 
instant. In the sequel we shall simply write W^ os (q) instead of W^ os (|_<zJ) when q is not an 
integer. If customer — m is still in service at time 0, we may write 

W^ os = V- m + B 

—m T— m + W^ os (A(T- m ,T- m + V- m + B- m )), 
where A(s, t) denotes the number of arrivals between times s and t. We obviously have 
W^ros — B—m — T- m + W^Q^(A(T- m , T- m + B- m )), 

+ V- m + B- m )). (C.l) 

The first bound gives 

oo 

P (Wros > x > B -m > (ma + x)(l - p), V- m < T- m < B-m + V- m < T- m + a;) 

m=l 

oo 

> E P ( B - m ~ T - m + W ROs(MT-m, T- m + B-m)) > x, B- m > (ma + x)(l - p) , 

m=l 

V—m ^ T— m < B—m, ~\~ V—m ^ T— m ~\~ x) 

oo 



> (1 - 5f V P(B- m - ma(l + e) + W^ os (-(l - e)B- m ) > x, B- m > (ma + x)(l-p), 
* — * a 

m=M 

1+e) < B-m < ma(l - e) -K + x), (C.2) 



ma 
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where, for fixed e > 0, 5 > 0, we have chosen K > such that P(V > K) < 5 and M > 
such that Ma(l — e) > K and, for all m > M and y > Ma(l — p), 

P ^(1 - e)ma < T_ m < (1 + e)ma, A(T. m , T_ m + y) > * y(l - e)^ > 1 - 6, 

which is possible by the strong law of large numbers. 

Note that the summation in (C.2) is actually truncated at m = (x — K)/(2e). For notation 
it is convenient to make the summation run from m = to oo. Adding the terms for m < M 
in (C.2) causes an error of the order o(P(B^ w > x)) and since 



^2?{ma{l-e) + x-K <B^ m <ma + x) < eO(P(B fw > x)), 

m=0 

oo 

P(B- m > {ma + x){l- p),ma < B_ m <ma(l + e)) < eO{P(B fw > x)), 



m=0 

we have, 



Yl P Wos > x, B_ m > (ma + x)(l-p), V. 

m ^ T— m < B— m -\- V— m — T— m + Xj 



m=l 

+ 7T^V2°( P (B fw > *)) 



(l-5f 

> V P(5_ m - ma(l + e) + Wros(-(1 - e)£-m) > s, £_ m > (ma + x)(l - p), 



m=0 



ma < -B_ m < ma + x 



Furthermore, by Lemma 3.1 we can find xq such that for all x > xq: 

rma+x / i \ 

/ dP{B < z) P I Wro S (-(1 - > x - z + ma(l + e) ) 



i 



z=max{(ma+3;)(l-p).mtt} V V P(l £j£ 



=max{(ma+x)(l— p),ma} \ V 

where 7 > depends one. In the last step we used that, as e — > 0, 

(1 — p){x — z + ma(l + e)) (1 — p)(x — z + ma) 
p(l — e)z pz 

uniformly in z within the area of integration. (This can be seen, using that z > (ma+x)(l— p) 
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and z > ma.) So that, with 5 — > 0, e — > and then 7 — > we have: 

00 

P (^ROS > x, B_ m > {ma + x)(l - p), V- m < T_ m < B- rn + V^ rn < T_ m + x) 



m=l 

m= g z=max{(ma+:r)(l — p),mo} V / 

+o(P(£ /w > x)). 

Next we derive an upper bound, 

00 

Yj P (^ROS > #-m > (ma + x)(l - p), < T_ m < £_ m + y„ m < T_ m + x) 

m=l 

00 

= p + " T - m + w Bos(MT- m , r_ m + + B_ m )) > x, 

m=l 

£_ m > (ma + x)(l - p), y_ m < T_ m < 5_ m + V- m < T_ m + x) 

< P{V > K)0(P(Bf w > x)) 

00 

+ Y P ( K + B -m ~ T_ m + W^o S (A(T_ m , T_ m + ^ + S_ m )) > x, 

m=l 

S_ m > (ma + x)(l - p), T_ m - if < 5_ m < T_ m + x) 

< (5 + P(F > if)) 0{P{B fw > x)) 

,1 +£, 

a 



+ V P(K + B_ m - ma(l - e) + <os(— (^ + £-m)) > x, 

-B-m > {ma + x)(l — p),ma{\ — e) — K < B- m < ma{l + e) + x^j , 
where, for fixed e > 0, we have chosen M > such that, for all m > M and y > Ma(l — p), 

P ((1 - e)ma < T_ m < (1 + e)ma, A{T_ m , T_ m + y) < ^y{l + e)^j > 1 - S. 

As in the lower bound we may let the summation run from m = 1 to 00 and replace the 
condition ma{l — e) — K < B- rn < ma{l + e) + x with ma — K < B- m < ma + x — K; 
the error we make is of the order eO{P{B^ w > x)). Also, replacing B > (x + ma){\ — p) by 
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B > (x + ma) (I — p) — K does not decrease the probability. 

oo 

P (^ROS > x, B„ m > (ma + s)(l - p), V- m < T_ m < B„ m + V- m < T_ m + x) 
< (e + S + P(V > K)) 0(P(B fw > x)) 

OO j 

+ Y P ( K - ma ^ - £ ) + B — + ^ROsl — (K + B-m)) > 



a:, 

a 

m=0 

B_ m > (ma + x)(l — p) — K, ma — K < B- m < ma + x — K 
= (e + 5 + P(V > K)) 0(P(B fw > x)) 

+ / dP(£ < z - K) p(Wr OS (—^z) >x + ma(l - e) - z) 

m _0 •/ z=max{(ma+i)(l- p),ma} \ Q 

< (e + 5 + P(V > K)) 0(P(B fw > x)) 

+ d + 7)Er +i dP (fl <,- J offi- (i - rt(,+mg -' ) )l 

m _Q J 2=max{(ma+x)(l— p),ma} y\ / / 

In the last step we use Lemma |3.1| and the uniform convergence of 

(1 — p)(x + ma(l — e) — z) 
p{l + e)z 

as e — > (7 > depends on e). 
Using that 

00 

P(max{ (ma + x)(l — p), ma} — K < B < max{ (ma + x)(l — p),ma}) = o(P(B^ w > x)), 

and 

00 

P(ma + x- K<5<ma + x) = o(P( J B /u ' > x)), 

m=0 

we may replace dP(-B > z — K) by dP(-B > z). Now let IT — > 00, e — )■ 0, 5 — >• and then 
7 — >• to conclude that 

00 

Y P (^ROS > x , B -m > {ma + x)(l- p), V- rn < T_ m < 5_ m + V^ m < T_ m + x) 

m=l 

< o(P(^- >,)> + £ r +I dP (B < „ (1 - (wx- +"«-»>) r 

m _Q ^ z=max{(ma+i)(l- p),ma} V / 

Together with the lower bound, this shows 

00 

Y P ( W *OS > x > B -m > (ma + x)(l-p), V_ 
„(P( B '» > ,)> + £ P + * dP(B < ,) (l - (1 -«)(- + "■"-»)) r 

m— o z=max{(mo+x)(l— p),ma} V P^ / 
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m=l 



The second and third term in Theorem 3.2 now readily follow, using that 



171 = 



q i 2=max{(m«+i)(l-/)) 1 mtt} V P^ 



v=0 J z=ma,x{(va+x)(l—p),va} V 

+ o(P(S- fu ' > x)). 

Part III. Finally, we deal with the last possible scenario in which service of the large 
customer ends before time 0. Thus, V- m + B_ m < T_ m . Suppose that customer — m + N 
is in service at time 0; 1 < N < m — 1. We can bound the waiting time of customer from 
below by 

Wrqs > Wi os (A(-T^ m ,0)-N) = W^ os (m-N), 

and from above by 

W^os < B^ m+N + W^ os (m-N + l). 

We start with the lower bound. We shall denote the number of departures in the interval 
[u,v) by D(u,v). Note that N = D(-T- m + V- m + B- m ,0) < D(-T- m + B_ m ,0). In the 
following we take e, 5, M and K such that P(V > K) < 5 and for all m > M, y > K, 

P ^(1 - e)ma < T_ m < (1 + e)ma, D(-y, 0) > (1 - e)^j > 1 - 5. 

We have, 

oo 

P ( W ROS > x,B-m > (ma + x)(l - p),V- m + B„ m < T_ m ) 

m=0 

oo j 

> (i_ 5 )2 £ p(V^ os (m -((l + £ )ma-5_ m ))>x ) 



m=M 



B- rri > (ma + x)(l — p), K + B- m < (1 — e)ma 

oo 1 

= eO(P(B*° > x)) + (1 - 5) 2 P(w^o S (m - (1 - e) - ((1 + e)ma - B_ m )) > x, 

m=0 ' 

B- rn > (ma + x)(l — p), B_ m < ma^j 

= eO(P(B fw > x)) 

OO /* TYl, Q. 1 

+ (l-5) 2 Y dP(B<z)p(w^ os (m-(l-e)-((l + e)ma-z))>x) 

m=0 Jz=(ma+x)(l-p) v P ' 

> eO(P(B^ w > x)) 



+ (1 - 7)(1 - *) 2 E f ° dP(5 < z) ( 1 - ) 

Jz=( ma +x)(i-p) v ^ - - p) J 



i 
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where 7 > depends on e. In the last step we used Lemma 3.1 and the uniform convergence 
of 

(1 ~ P)x 

/3(m - (1 - e)i ((1 + e)ma - z)) 

as e — > 0. Similar to Part II it can be shown that 



tr dP(B<;)(l- "-f - Y* 

^ Q Jz=(ma+x)(l- P ) \ z-ma{\- p)J 

= (p(*/« > x)) + r r &?(B <z)(i- { \-j> x ) 

Jv=oJz=(va+x)(l-p) V Z-VOt(l-p)J 



1 



I v=0 J z={va+x){l— p) 

Letting e — > 0, <5 — > and 7 — ?■ we thus have proved that 



2 p (Wros > s -m > (ma + x)(l- p),V_ m + 5_ m < T_ m ) 

m=0 

> (p(b*° > x)) + r r° dP( J B < Z ) (1 - {i ~ p ^ x ) 

J v =0 Jz=(va+x)(l-p) V z ~ - PJ / 



1 

l-p 



>v=0 J z=(va+x)(l— p) 

It remains to show that the right-hand side is also an upper bound for the left-hand side. 
Recall that N = JD(-T_ m + V- m + B_ m , 0) and W^ os < B_ m+N + W^ os (m - N + 1). Note 
that, if e > and J > and M such that P(TL m > (1 + £)ma) < 5 for all m > M, then 

00 

P (-B_ m > (ma +x)(l - p), V. m + £_ m < T_ m , V. m + 5_ m > (1 - 2e)ma) 
< (5 + P(V > K))0{P(Bf w > x)) 

+ ^2 P ( B -m > (ma + x)(l - p), (1 - 2£)ma - K < B_ m < (1 + e)ma) 

= (e + 5 + P(V > K))0(P(Bf w > x)). 

We shall use this in what follows. In addition, let M and K be such that P(V > X) < 5, and 
for all m > M and 7/ > eMa, 

P(£>(-y,O)<(l-e)i//0)<5. 
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Also, we take L such that P{B > L) < 5. 

oo 

E P ( W bos > x, £- w > (ma + x)(l - p), V- m + £- m < T- m ) 

m=0 

< (e + 5 + P(V > K) + P(fl > L))0(P( J B /w > s)) 

oo 

+ E P (^ROsC™ - JV + 1) > x - L, B^ m > (ma + x - L)(l - p), V- m + £_ m < (1 - 2e)ma) 

m=M 

< (e + 8)0(P(Bf w >x + L)) 

OO j 

+ E P(w / ROs(^-(l-£)^((l-^«-^-^-m) + l)>X, 
m=M ^ 

S-m > ( 7Tia + ~~ p),B^ m < (1 — 2e)ma^ 
= (e + 5)0(P(B*°>a;)) 

+ E P (<os(™ " (1 " " f)« " ^ " B_ m ) + 1) > x, 

m=0 " 

B-m > (ma + x)(l — p) — K — 1, B_ m < ma — K — 1^ 
= (e + 8)0(P(Bf« > x)) + (1 + 7 ) /°° f " dP(B <z)(l- {l ~ ^ . ) . 

As before, in the last step we use Lemma |3.1[ the uniform convergence of 

(1 ~ P)x 

/3m - (1 - e)((l - e)ma -K - B_ m ) + 0' 

as e — > and the fact that replacing the summation with an integral and dP(i? < z — K — 1) 
with dP(B < z) introduces an error of the order o(P(B^ w > x)). Letting e — > 0, 6 — > and 



7—7-0 yields the last term in Theorem 3.2 □ 



D Random and deterministic arrivals 

The following lemma states that when interested in events involving a large service time, we 
may in fact ignore the randomness in the arrival process and replace it by a deterministic 
arrival process with the same mean arrival rate. Thus, heuristically, we may concentrate on 
the D/G/l queue instead of the GI/G/1 queue. Although the lemma is not explicitly used 
in the paper, it has been very useful in guiding us to the proof of several of its results. We 
formulate it here because we expect that a reduction to a deterministic arrival process will 
often be helpful in proving tail asymptotics, see also Baccelli and Foss [2]. 

Lemma D.l. If B = Bo has a finite first moment and if its integrated tail distribution belongs 
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to C, then, for any constant c > 0, as x — >• oo, 

oo oo 

P {B > x + cT_ m ) ~ ^ P(B > x + cT_ m ,B > x + 



oo 

" " " T ' ' cma) 

m=0 m=0 



]T P(£ > x + 



oo 

cma) 

m=0 

~ ( />' > ./• i . 



Proof. From Appendix |A| (Property (lc) in Appendix |Aj), the integrated tail distribution of 
f? belongs to £ if and only if B^ w G £. 

Lower bound. For any e£ (0, 1), choose i? > such that 

inf P (T_ m < ma(l + e) + R) > 1 - e. 

m>0 

Then, as x — > oo, 

oo oo 

J2 p ( B > x + cT-m) > J2 P ( B>X + cma ( l + e ) + cR i T -m < ma(l + e) + R) 

m=0 

oo 

> (1 -e) P ( B > x + cma{\ + e) + cR) 

> ^— -£p(V w >x + cR 
1 + e c V 

1_ ^pfi?^>x|. 



m=0 m=0 

oo 



1 + e c 

Letting e \. 0, we get the right lower bound. 

Upper bound. Fix any e E (0, 1). Since T_ m are partial sums of non-negative i.i.d. r.v.'s 
with a finite positive mean a, 



K = ^ P (T- m < ma(l - e)) < oo. 

Then, as x — > oo, 



m=0 



^ P (B > x + cT_ m ) < ^ P (B > i + cm«(l - e),T_ ra > ma(l - e)) 

oo 

+ ^ P(B > x,T_ m < ma(l - e)) 

oo 

< J2P(B>x + cma(l-e)) + KP(B>x) 



m=0 m=0 

oo 



m=0 



< — — ^P (B fw > x-ca) + KP(B > x) 

1 — £ C V / 

1 P p(b^>x 



1 — £ C 

since S^ 1 " G C and P(5 > x) = o(P(B fw > x)), see Appendix [X] (Properties (lb) and ([lc])). 
Letting e 1 0, we get an upper bound which coincides with the lower bound. 

□ 
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